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Abstract 

The generalized Tower of Hanoi problem with h > 4 pegs is known to require a sub-exponentially 
fast growing number of moves in order to transfer a pile of n disks from one peg to another. In this 
paper we study the Path^ variant, where the pegs are placed along a line, and disks can be moved 
from a peg to its nearest neighbor(s) only. 

Whereas in the simple variant there are h{h—l)/2 possible bi-directional interconnections among 
pegs, here there are only /t — 1 of them. Despite the significant reduction in the number of inter- 
connections, the number of moves needed to transfer a pile of n disks between any two pegs also 
grows sub-exponentially as a function of n. 

We study these graphs, identify sets of mutually recursive tasks, and obtain a relatively tight 
upper bound for the number of moves, depending on h,n and the source and destination pegs. 

Keywords: Tower of Hanoi, path graphs, analysis of algorithms 

1 Introduction 

In the well-known Tower of Hanoi problem, proposed over a hundred years ago by 
Lucas [20], a player is given 3 pegs and a certain number n of disks of distinct sizes, 
and is required to transfer them from one peg to another. Initially all disks are stacked 
(composing a tower) on the first peg (the source) ordered monotonically by size, with the 
smallest at the top and the largest at the bottom. The goal is to transfer them to the 
third peg (the destination), moving only topmost disks, and never placing a disk on top 
of a smaller one. The well-known recursive algorithm that accomplishes this task requires 
2" — 1 steps, and is the unique optimal algorithm for the problem. The educational 
aspects of the Tower of Hanoi puzzle have been reinforced recently, by a series of papers 
1 Research supported in part by the Sapir Academic College, Israel. 



by Minsker ([231 [211 125]), composing variants for the sake of studying their combinatorial 
as well as algorithmic aspects. 

Work on this problem still goes on, studying properties of solution instances, as well 
as variants of the original problem. Connections between Pascal's triangle, the Sierpihski 
gasket and the Tower of Hanoi are established in [16], and to some classical numbers 
in [I^. In [I] it is shown that, with a certain way of coding the moves, a string which 
represents an optimal solution is square-free. This line is extended in [2J. Another 
direction was concerned with various generalizations, such as having any initial and final 
configurations [T4], assigning colors to disks (cf. [21] and [22j for recent papers on the 
subject), and relaxing the placement rule of disks by allowing a disk to be placed on top 
of a smaller one under prescribed conditions [^ ITUl [TT] . 

A natural extension of the original problem is obtained by adding pegs. One of the 
earliest versions is "The Reve's Puzzle" [121 PP- 1-2]. There it was presented in a 
limited form: 4 pegs and specified numbers of disks. The general setup of the problem, 
with any number /i > 3 of pegs and any number of disks, was suggested in [2Bj, with 
solutions in [2^ and [I^, shown recently to be identical [TT]. An analysis of the algorithm 
reveals, somewhat surprisingly, that the solution grows sub-exponentially, at the rate of 
Q{^/n2'^^) for /i = 4 (cf. [30J). The lower bound issue was considered in [32J and [5], 
where it has been shown that the minimal number of moves grows roughly at the same 
rate. 

An imposition of movement restrictions among pegs generates many variants, and 
calls for representing variants by digraphs, where a vertex designates a peg, and an 
arc represents the permission to move a disk in the appropriate direction. In [31 [13], 
the uni-directional cyclic 3-peg variant (Cyclica) has been studied, and the average 
distance between the nodes - in [21]. In [27], the "three-in-a-row" arrangement (Paths) is 
discussed. A unified treatment of all 3-peg variants is given in [2S]- The (uni-directional) 
Cyclic4 is discussed for the first time in [2Z], and [20] studies other 4-peg variants: 
Star4 and Path4, presenting a sub-exponential algorithm for Star4. The Cyclic/j for 
any number of pegs h> A has been studied in [5] and proved to be exponential for any 
specified h. Identification of the longest task, for certain variants, has been resolved in 
[7J. 

The only requirement for the problem to be solved for any number of disks is that 
the variant is represented by a strongly-connected directed graph. An interesting line of 
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research has been taken in [12], [1], and |5], where non-strongly-connected graphs are 
being studied. 

In this paper we study the Path/^ variant, with a fixed number /i > 4 of pegs, whose 
complexity issue has been left open. We devise an efficient algorithm which moves a 
column of n disks between any pair of pegs, and supply an explicit subexponential upper 
bound on the number of moves, for each h. 

Notations and definitions are given in Section [21 the main results in Section [31 the 
proof of the 4-peg case in Section [4] and that of the general case in Section [51 

2 Preliminaries 

We study the Path/j (a.k.a. h-in-a-row) variant, with a fixed number /i > 4 of pegs. 
We denote the pegs of Path/^, from left to right, by 1, . . . , /i. Let the sizes of the disks 
be 1,2, ... ,n. For convenience, we identify the name of a disk with its size. 

For the statements and algorithms of the paper, it is required to introduce the notion of 
a block — a set of disks of consecutive sizes. The minimum (respectively, maximum) size 
of a disk in a block B is denoted by i?min (resp., i?max), and the number i?max — -Bmin + 1 
of disks in i? — by \B\. A block B is lighter than another block B' if -Bmax < -Bmin- 

A configuration is a legal distribution of the n disks among the h pegs. A perfect 
configuration is one in which all the disks reside on the same peg. Such a configuration 
is denoted by Rh.i.n, where h is the number of pegs, i the peg holding the disks, and n 
the number of disks. 

For a sequence of moves M, henceforth move-sequence, we denote by the reverse 
move-sequence, comprising the moves that cause the reverse effect. That is, the order 
of the moves is reversed and each move of the original sequence is reversed. Clearly, if 
applying M to configuration Ci results in reaching configuration C2, then applying 
to C2 results in configuration Ci. (Note that this is true if and only if the peg structure 
is a graph; for digraphs in general this is not true.) 

A problem instance, henceforth a task, is given by a pair of configurations, an initial 
configuration Ci and a final configuration C2, where we are required to move from Ci 
to C2 in a minimal number of moves. The task, as well as a minimal-length solution of 
it, is denoted by Ci C2, and the minimum number of moves needed to get from Ci 
to C2 is denoted by |Ci — )► C2I. 
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In this paper we focus on perfect tasks — problem instances whose initial and final 
configurations are both perfect. The peg associated with the initial (respectively, final) 
configuration of a perfect task is naturally referred to as the source (resp., destination). 
Clearly, for any positive integers h,n and 1 < i < j < h, we have \Rh,i,n — ^ Rh,j,n\ = 
\Rh,j,n — > Rh,i,n\- We shall henceforth restrict our attention in the sequel to tasks in 
which the source peg is to the left (i.e. has a lower peg index) of the destination peg. 

For n > 1, denote by Path(/i,n) the minimal number of moves which suffices for 
transferring a block of size n between all pairs of perfect configurations in Path^j, namely, 

Path(/i,n) = max \Rh,i,n Rh,j,n\ ■ 

l<i<j<h 

For a real number x, let round(x) be the integer closest to x (where round(x) = [x] 
for X = n + 0.5). For a pair of positive integers p and q, with p < q, we denote the set 
{p, . . . ,q} by [p, q], and [1, . . . , by [q]. In what follows, we do not distinguish between a 
move-sequence and an algorithm generating it, if this does not lead to a misunderstanding. 

3 Main results 

The main question the paper addresses is: what is the complexity of Path(/i,n)? 
An upper bound is provided by 

1 

Theorem 3.1 Path(/i, n) < Chu"'^ ■ 3^''"''"' , for allh>3 and n, where: 

'^h = - 

Oh ' y ~ 32-(l/30)l/3 )■ 

In particular, Path(/i,r?,) grows subexponentially as a function of n for h>A. 

Of course, as a lower bound for Path(/i, n) one may use any lower bound for the 
number of moves required to move a tower of size n from one peg to another over the 
complete graph on h vertices. By [S], such a lower bound is given by 2(^+°(^))("(^~2)') ''"^ , 
which is "not very far" from our upper bound for Path(/i,n). 

The following theorem identifies the hardest perfect task for the particular case h = A. 
It also provides a tighter upper bound for Path4 than the one given in Theorem 13.11 

Theorem 3.2 For every n > 1; 
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(a) \R4,i,n ^4,j,n| < |i?4,l,n. ^ ^4,4,n,| foT 1 < i < j < 4:,{i,j) ^ (1,4). In 

particular, Path(4, n) = |-R4,i,„ — )■ -R4,4,n|- 

(b) Path(4,n) < l.e^nS^. 

4 Proof of Theorem 13.2 

4.1 On the relation between various tasks in Path4 

We start with a result of some independent interest, which holds for general h. 

Lemma 4.1 Let C be a configuration with n > 1 disks, arranged arbitrarily on pegs 
1, . . . , h — 2, with pegs h — 1 and h empty. Then: 

,h—2,n\ < \C- y Rh,h,n\i {C ^ Rh,h-l,n\ < \C ^ Rh,h,n\ ■ 

Proof: We detail the proof for \C — )■ Rh,h-2,n\ < |C — Rh,h,n\- The proof for the 
second inequality is similar. 

The proof is by induction on n. The basis n = 1 is trivial. Let n > 2, assume that the 
statement holds for up to n — 1 disks, and let C be a configuration as in the statement 
of the lemma and M a move-sequence transferring from C to Rh,h,n- Before the last 
move of disk n (to peg h), a configuration C , in which all n — 1 disks 1, 2, . . . , n — 1 
are distributed among pegs 1, . . . , h — 2, is reached. Let M' (respectively, M") be the 
subsequence of M, consisting of all moves that come before (resp., after) the last move of 
disk n. Notice that M" transfers from C (considered as a configuration of n — 1 disks) to 
Rh,h,n-i- By the induction hypothesis, there exists a move-sequence M^'_2 that transfers 
from C to Rh,h~2,n-i, which is strictly shorter than M" . Let M^_2 be the move-sequence 
obtained from M' by omitting all moves disk n makes after reaching peg h — 2 for the 
first time. Concatenating M^_2 with M'^_2, we obtain a legal move-sequence, strictly 
shorter than M, transferring from C to Rh,h-2,n- The required result follows. 

Due to symmetries, there are actually only four essentially distinct perfect tasks in 

Path4: i?4^i,„ ~^ RA,2,n, RA,l,n -^4,3,™. RA,l,n ~> RA,A,n: 3rid i?4,2,n RA,3,n- By 

Lemma 14711 taking h = A and C to be various perfect configurations, we obtain for any 

n>l 

• 1-^4,1,71 — ^ RA,2,n\ < 1-^4,1,71 — ^ -^4,4,71 1- 

• 1-^4,2,71 — ^ RA,3,n\ < 1-^4,2,71 — ^ -R4,4,n| = |-R4,l,n -R4,3,n| < 1-^4,1,™ — ^ RA,A,n\- 
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These inequalities establish part (a) of Theorem 13.21 

In Table 1 we present the (distinct) numbers iR^^i^n ^4,j,n| for 1 < n < 11. The 
entries have been calculated by finding the distance between the vertices i?4,i,„ and i?4j> 
in the graph of all configurations of n disks on Path4 using breadth-first search. 
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1 ^ 3 


1 ^ 4 
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0.634 
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10 


0.786 
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7 


9 


12 


19 


0.744 


4 


14 


18 


22 


34 


0.760 


5 


23 


29 


36 


57 


0.790 


6 


34 


44 


54 


88 


0.799 


7 


53 


69 


78 


123 


0.762 


8 


78 


96 


112 


176 


0.768 


9 


105 


133 


158 


253 


0.798 


10 


138 


182 


212 


342 


0.795 


11 


187 


241 


272 


449 


0.783 



Table 1: The minimal numbers of moves for the 4 different perfect tasks in Path4 
The table prompts 

Question 1 Is it the case that |i?4,i,ri — ^ -R4,2,n| < |-R4,i,n — ^ -R4,3,n| and |-R4,2,n 

RA,3,n\ < \R4,l,n RA,2,n\ foT all n > 3? 

Both of these inequalities seem intuitively quite plausible. 
4.2 Upper bound for Path(4, n) 

In this subsection we present the algorithm FourMove for moving a block B of size n from 
peg 1 to peg 4 in Path4, requiring no more than 1. 6^713^^ moves. By Theorem 13.2( a). 
this will imply Theorem 13.2( b). The description of FourMove is given in Algorithm [T] 

Prior to its main stages, FourMove partitions B into three: a block containing the 
smallest disks, a block containing the larger ones, and a block containing a single disk 
- the largest one. These blocks are denoted Bg, Bi, {Bj^ax} respectively, with m = 

|-BzU{^max}|- Thus 5, = [Smm,5max- W-] and = [^max - "I + 1 , S^ax " 1] ■ In the 

three principal stages that follow: Spread, Circular shift and Accumulate, the moves are 
done based on these blocks. In Spread, Bg is transferred to the farthest peg - number 4, 
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Algorithm 1 FourMove(i?) 



/* _B is a block of disks, partitioned into BgU BiU {i?max}, where Bg is the set of smallest */ 
/* disks, Bi — the larger, {-Bmax} — the largest. ThreeMove(i3,s,d,a) returns the shortest */ 
/* move-sequence for moving B from s to d using a, referring to the graph induced by these */ 
/* three vertices as Paths. A move of disk n from peg i to peg i + 1 is denoted by */ 
/* The '*' denotes concatenation of move-sequences. */ 

M ^ [] /* an empty sequence */ 

if B 7^ then 
n ^ \B\ 

m -s— round(\/2n) 

Bg ^ [Bnun, B^i^y^ — to] 

Bl ^ [-Bmax - TO + 1, Bmax " 1] 

Ms FourMove(Bs) 
/* Spread: */ 

M ^ Ms* ThreeMove(B/, 1,3,2) * ti^2,n 
/* Circular shift: */ 

M M *M-^* ThreeMove(B,,3,4, 2) * t2,3,„ * ThreeMove(B/, 4, 2, 3) 
/* Accumulate: */ 

M ^ M* t3,4^„ * ThreeMove(B,, 2, 4, 3) * Ms 
end if 
return M 



Bl to peg 3, {Smax} to peg 2. In Accumulate, the opposite is done: these blocks are 
gathered on the destination peg. In-between the algorithm performs the Circular shift 
stage, whose role is to reverse the order of the blocks, so that it will be possible to 
perform the Accumulate stage. It is easy to verify that, as the execution of the algorithm 
terminates, all the blocks are legally gathered on the destination peg 4, as required. 

The ThreeMove procedure in Algorithm [T] produces move-sequences for Bi using only 
three (consecutive) pegs, which is exactly as moving it in Paths. To this end, we use the 
algorithm of [27j, which transfers a block in a minimal number of moves between any 
two pegs in Paths, requiring 3" — 1 moves to transfer n disks between the two farthest 
pegs, and half that number of moves between neighboring pegs. 

Denote by T(n) the number of moves required by FourMove for a block of size n, and 
define T(0) = 0. Each of the three recursive invocations of the algorithm FourMove with 
Bs requires T{n — m) moves. Observe that, for a positive integer n, we have 1 < m = 
round(A/2n) < n. Employing the abovementioned results regarding the number of moves 
required by ThreeMove, it is easy to see that the total number of moves required by Bg 
is 3(3™-^ - 1) + 1(3'"-^ - 1). Finally, {S^ax} performs 3 moves. Altogether, for n > 1 
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we have: 



T(n) 



3 ■ T(r2 - m) + I ■ (3 



m— 1 



1) + 3 



3 • r(ri - m) + I ■ 3™ 



(1) 



where m = round(\/2n). 

Next, we prove by induction that T{n) < 1.6y/n- 3^^, implying the required result. 
For the induction basis, we note that the inequality has been verified manually for all 
values of 72 < 8. Let n > 9, and assume that the inequality holds when n is replaced by 
a smaller integer. To prove it for n, denote first /3 = m— y/2n. (Clearly, — | < /3 < |.) 
Note that 



1-^ < 1- 
n 



2/3 + 1 
4n ' 



(2) 



which can be verified by squaring both sides of the inequality, noting that the right-hand 
side is positive for n > 2. Thus, by the induction hypothesis and (EI), 

T{n) = 3-T{n-m) + l-3"' -I 

D 2 



< 3■1.6^/■ 



n — m 



6 



3-1.6v^,/l---3^v^+--3"^ 
n 6 



< 3-1.6v^ 1 



}x/2n 



/2n 



1 



2n 



2/3 + l^3v/2;^(l--^-M±i; 



An 



i.eVnJ^i- 

1.6v/n(l- 



1 


2/3 + 1 




An 


1 


2^ + 1 


\/2n 


An 



2/3 + 1 7 a 

3 2v^ + _ . 

6 



6'^ 



e 2^ H — ■3'^ 



7 
6 



Notice that the term — |^^ln3 is negative for /3 > 



Hence, using the fact that 



<1 + X + for a; < 0, we obtain 
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Tjn) 
1.6-3^ 



^ AT/'l _ 1 2/3+A 2/3+1 . o I (2/3+l)^ln^3 N , 35 o0 



(2^+1) In 3 
2v^ 



I / 2/3+1 I (2/3+1) In 3 , (2/3+1)^ In^ 3 

, / (2/3+l)^In3 (2/3+l)^ln^3 (2/3+1)^ In^ 3 \ , 35 o/3 
8V2-n 16v/2-n 64nx/H J 48 ' 

(vs-A-^) 

I ( 2^+1 , (2^+1) in 3 , (2/3+1)^ In^ 3 , (2/3+1)^ In^ 3 \ 
\^ 4v^ ~^ 16xAI 16\/2-n / 

, / (2/3+1)^ In 3 _ (2/3+1)2 In^ 3 _ (2/;+l)3 In^ 3 \ , 35 o/3 
~^ 8V2-n 8v^-n 64n0^ J 48 ' " 



Observe that 

(2/3+ 1)^ In 3 _ (2^ + l)^ln^3 _ (2/3 + l)3ln^3 
8\/2 • n 8^2 • n 64 • n^ri 

for ^ > — |, and therefore 



< 



Tjn) 
1.6-3^ 



I / 2/3+1 I (2/3+1) ln3 . (2/3+1)^ In" 3 . (2/3+1)^ In^ 3 \ , 35 q/3 
4VH 16V2-n ^ 48 ' 

= (Vi-;J5-«SF) 

+ fUs (4(ln3 - 1) + (2/) + 1) In^ 3(1 + ^))+!. 3". 
Now for /3 > — I the expression 

+ ^ 4(ln3 - 1) + {2/3 + 1) W 3 f 1 + (3) 



increases as a function of /3 and decreases as a function of n. Hence its maximal value 
in the range — | < /3 < |, n > 9, is obtained for /3 = \,n = 9. Since its value at that 
point is less than 0.15, we have 
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1.6 -Sv^ V 72 2V2 48 y ^ ' 

It is easy to verify that for — | < /3 < |, 

is maximized at /3 = |, and /(|) < 0. Consequently, the right-hand side of (jH) is smaller 
than ^/n, and we are done. 

4.2.1 A better upper bound for Path(4, n) 

We reduced the problem of upper bounding Path(4, n) to the problem of upper bounding 
the following recurrence formula, which might be of independent interest: 

/ 7 1 

T{n) = mill 3 • r(n - m) + - ■ 3™ - - 

l<m<n \ 6 2 



We obtained an upper bound of 1.6 ■ y/n ■ 3^^" for this recurrence formula, which is tight 
up to the leading constant 1.6. 

One way to decrease the constant 1.6 is to show, using a computer program, that a 
better upper bound holds for all values of n < n' , for some huge integer n'. This will 
serve as a significantly more elaborate induction basis than the one that we use above 
(i.e., n < 8), and consequently, it would suffice to prove the induction step for n > n' 
only. The maximum value of the function g(n, (3) defined in (ED for a huge integer n > n' 
and — I < /3 < I is some tiny number e = e{n), and so, substituting 0.15 with e in ([5]) 
yields: 

with /(/?) - /'(/?) = 1.6(0.15 - 5) ^ 1.6 ■ 0.15. The difference 1.6(0.15 - e), between 
f{f3) and /'(/?), for a sufficiently small e, enables us to decrease the leading constant 1.6 
to approach 1.365, yielding an upper bound that approaches 1.365 ■ y/n ■ 3^^. 

A more involved method for decreasing the above constant is to choose another value 
for m. For technical convenience, we fixed m = round(-\/2n), but this choice of m is 
inherently suboptimal. By following the method outlined in the previous paragraph, and 
setting m = \_\/2n + (l — i^)J. one can achieve a constant that approaches 1.105, 
yielding an upper bound of approximately 1.105 ■ y/n ■ 3^^. 



10 



5 Proof of Theorem 13.1 



The proof of Theorem l3.ll is organized as follows. Generally, we would like to show how 
one can move a column of n disks from any source peg s to any destination peg d such 
that the number of moves is bounded above as the theorem states. For simplicity, we 
start by presenting an algorithm for the case where s = l,d = h. This will be done in 
Section [5H Then we present an algorithm for the general case (Section [5.20 . We note 



that, in fact, the first algorithm does employ the second. An important point in both 
cases is a partitioning of the set of disks to blocks, which will be discussed in Section 1531 
Time analysis of the two algorithms will be provided in Section [5341 

5.1 Moving disks between the farthermost pegs 

Here we present FarthestMove (Algorithm [21), designed to move a block B of n disks 
between the two farthest pegs in Path/j, where h > 3. 

We partition B in some way to blocks Bi(h, B), B2{h, B), . . . , Bh-i{h, B) of disks. 
Whenever h and B are implied by the context, we write Bi instead of Bi(h,B). The 
block Bi consists of the smallest ni disks 1,2, .. . ,ni, the block B2 — of the n2 next 
smallest disks ni + 1, ni + 2, . . . , ni + ri2, and so forth. Similarly to the shorthand used 
when denoting blocks, we may write Ui (with a possible superscript) instead of 
For any i e [h - 1], let B{i) = [j^zl Bj and n{i) = \B{i)\ = Sj^/n^, where uj = \Bj\. 
(Note that B{1) = B and n(l) = n.) 

The determination of the sizes rij is crucial for the number of moves the algorithm 
makes, and will be explained later. However, for the algorithm to work correctly, it is 
only required for Bh_i to consist of the single disk Smax — the largest. The algorithm 
consists of three phases (see Figure [T]for an illustration): 

• Spread: Move the h — 2 {= d — s — 1) first blocks Bi, . . . , Bh_2 fi'orn the source 
peg s to pegs d, . . . , d — h + 3, respectively. It consists of h — 2 iterations. At the 
j-th iteration, j E [h — 2], block Bj is (recursively) moved from s to d — j + 1, using 
the set [1, d — j + 1] of available pegs. (Note that the 1-disk block B^^i has not 
been moved from s to s + 1. It is more convenient for us to view this move as the 
first move of the next stage.) 

• Reverse: The role of this phase is to reverse the positions of the h — 1 blocks on 
the h pegs, i.e., a block residing, at the beginning of this phase, on peg s + i — 1 
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reaches, at the end of the phase, its reflected position — peg d — i + 1. The phase 
starts by moving the last block B^-i from s to s + 1. Then, h — 2 rounds are 
carried out, each of which brings the next larger block to its reflected position. The 
following highlights the way each round j achieves its goal: 

— Before this round, blocks Bi, . . . , Sj-i are on pegs s, . . . , s+j — 2, respectively; 
peg s+j — 1 is vacant; blocks Bj, . . . , Bh^i are on pegs d, . . . , s+j, respectively. 

— Block Bj is moved from to s + j — 1. 

— Blocks -Bj+i, . . . , -Bft-i are each shifted one peg to the right. 

— At the end of the round, blocks Bi, . . . , Bj are on pegs s, . . . , s + j — 1, respec- 
tively; peg s + j is vacant; blocks -Bj+i, . . . , -B/i-i are on pegs d, . . . , s + j + 1, 
respectively. 

Thus, as a result of this phase, block Bh-i is moved from s to d and, for each 
j E [h — 2], block Bj is moved from peg d — j + 1 to the reflected position, namely, 
peg s + i - 1. 

• Accumulate: The role of this phase is symmetrical to that of Spread, i.e., to 
move the h — 2 first blocks Bh-2, ■ ■ ■ , Bi from pegs d — 2, . . . , s, respectively, to d. 
Similarly, it consists of /i — 2 iterations, where at the j-th iteration block Bh-i-j is 
moved from s + /i — 2 — jtod using the set [s + /i — 2 — j, d] of available pegs. 

It is easy to verify that, as the execution of the algorithm terminates, all the blocks are 
legally gathered on d, as required. The formal description of the algorithm FarthestMove 
is given in Algorithm [2 

5.2 Moving disks between any pegs 

The general algorithm for moving a block of disks, between any two pegs s and d, in 
Path/j, is presented here. For convenience we assume that s < d. This does not effect 
the generality of the algorithm since, as was mentioned in Section [21 if M is a solution 
of Rh,s,n -> Rh,d,n, then M'^ is a solution of Rh,d,n Rh,s,n- 

The issue of partitioning the disk set is handled exactly as it was done in FarthestMove. 
Algorithm GeneralMove consists of five phases: two spread phases, a phase in which the 
remainder disks are moved, and two accumulate phases. The set of available pegs is 
denoted by A, and its smallest and largest pegs by A^am and Amax, respectively. 
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f-move 
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1 2 3 4 5 

f-move : A move of a block between the current farthest pegs. 

Done by a call to FarthestMove . 
a-move : An actual move, of the largest disk. 

shift : A move of a block between two neighboring pegs. 

Done by a call to GeneralMove. 

Figure 1: Main steps of FarthestMove for Paths, with B = Ui=i ^i' s — 1 and d= 5. 
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Algorithm 2 FarthestMove(i?, s, d) 



/* The procedure moves a block B from the leftmost peg s to the rightmost peg d. Prior */ 
/* to its main stages, the block is partitioned into h — 1 blocks which are treated as 'atomic */ 
/* units'. At the first of the main stages, these blocks are spread along the pegs; at the */ 
/* second - their order is reversed; at the third stage they are accumulated on the destination */ 
/* peg. The procedure requires that s < d. li d — s = 1, then \B\ < 1. The denotes */ 
/* concatenation of move-sequences. */ 

r [] /* initializing the result sequence */ 

it B then 
d- s + 1 

{Bi,...,Bh-i) = Partition(/i,B) /* Algorithm [5] below */ 

/* Spread: */ 

for j 1 to /i — 2 do 

/* At each step, the next block moves to the farthest available peg. */ 

T T * FartliestMove(Sj, s,d- j + 1) 
end for 

/* Reverse: */ 

T <— T * tB„,ax.s,s+i /* Moving the largest disk once a peg to the right. */ 
A/ <— [] /* Initializing the temporary move-sequence. */ 

for j •(— 1 to /i — 2 do 

/* Block Bj moves to the peg on which it will stay for the rest of this phase. */ 

Mf FarthestMove(i?j , s + j - l,d) 

M ^ M * MJ^ 

for i-i— j + lto/i — 2 do 

/* Each block whose index is higher than j is shifted one peg to the right. */ 
M ^ M * GeneralMove(Bi, d + j - i, d + j + 1 - i, [s + j, d + j + 1 - i]) 

end for 

T ^ T * M * if3„„„,j+i,j+2 /* Moving the largest disk a peg to the right. */ 
end for 

/* Accumulate: */ 

for j -f— 1 to ft. — 2 do 

/* At each step, the next block is gathered on the destination peg. */ 
T ^T* FartliestMove(Sft_i_j ,s + h-2- j,d) 
end for 
end if 
return T 
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• LeftSpread: In this phase the s — fifst blocks Bi, . . . , are taken from 
peg s to pegs Ainin, . . . , s — 1, respectively. It consists of s — A^in iterations. At 
the j-th iteration, I < j < s — A^i^, block Bj is (recursively) moved from s to 
^min + i - 1 using the set [Amm + i - 1, ^max] of available pegs. 

• RightSpread: Here, the Am^x — d next blocks are taken, from peg s to pegs 
^max, ■ ■ ■ ,d + 1, respectively. At each iteration j, where 1 < j < Amax — d, block 
^^-Amin+j '2 moved from s to Anax - i + 1, using [s, A^ax - j + !]■ Since at 
each iteration the source and destination are at the opposite ends of the currently 
available set of free pegs, the move is done using algorithm FarthestMove. 

• MoveRemainder: In this phase, the remaining B{d — s) blocks are moved from 
s to d. Since, as before, the source and destination are at the opposite sides of the 
set [s,d] of available pegs, this is done by algorithm FarthestMove. 

• LeftAccumulate: The role of this phase is symmetrical to that of RightSpread, 
that is, move B^-A^i^+i, • • • , Bs+\A\-i-d from A^ax, . . . , d + I to d, respectively. It 
consists of A„i^^ — d iterations, where at iteration j, block Bs+\A\~d-j is moved from 
d + j to d using [s,d + j]. Unlike RightSpread, the moves made in this phase are 
not between the two farthest available pegs. 

• RightAccumulate: This phase is symmetrical to LeftSpread, consisting of s — 
Amin iterations where, at the j-th iteration, Bs^j is moved from peg s — j to peg 

d, using [S - j, Amax]- 

It is easy to verify that, as the algorithm terminates, all the blocks are legally gathered on 
the destination peg d, as required (see Figure|2]for an illustration). The correctness proof 
is omitted. The formal description of GeneralMove is given in Algorithm [31 Note that, if 
the source and destination pegs are at the opposite sides of A, then GeneralMove does 
the same as FarthestMove. 

5.3 Partitioning the disks into blocks 

In this section we discuss how to set the sizes of the blocks such that the number of 
moves will be relatively low. The general idea is to view the h — 1 blocks as 'atomic' units, 
each occupying a single peg (except for when it is moved). During the process of moving a 
block Bi from one peg s to another peg d, the other blocks stay intact. Furthermore, the 



15 




RightSpread 



far-move 



MoveRemainder 




LeftAccumulate 



far-mov( 



^-^left-move 



B2 









RightAccumulate 



right-move 



-B4 



Figure 2: An illustration of the execution of the algorithm GeneralMove for Paths, with B = -Bj 
s = 2, d = A and A = [5]. 
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Algorithm 3 GeneralMove(B, s, d, A) 



/* The procedure moves a block B from any peg s to any other peg d. The partitioning is as in */ 
/* FarthestMove. Then some blocks are spread to the left of s and to the right of d; next - the */ 
/* remaining blocks are moved; eventually, the smallest disks are accumulated on the destination */ 
/* peg. The procedure requires that s,d & A, s < d. If d — s = 1, then \B\ < 1. The '*' denotes */ 
/* concatenation of move-sequences. */ 

r ^ [] /* an empty sequence */ 

if B ^ then 

= Partitioii(|A|,B) /* Algorithm [5] below */ 
/* LeftSpread: */ 

for j 1 to S — Aniin do 

/* At each step, the next block moves to the farthest available peg on the left side. Since */ 
/* the source is not necessarily at the rightmost available peg, this is done by GeneralMove.*/ 

fn ^ ^min + J ~ 1 

M ^ GeneralMove(Bj , /„, s, [/„, Amax]) 
end for 

/* RightSpread: */ 

for j ^ 1 to Aniax d do 

/* At each step, the next block moves to the farthest available peg on the right side. */ 
/* Since the source is at the leftmost available peg, this is done by FarthestMove.*/ 

fx ^ ^max ~ J + 1 

T ^ T ^ FarthestMove(_Bs__4^^.^ s, fx) 
end for 

/* MoveRemainder: */ 

T * FarthestMove(B(|^| -d + s), s, d) 

/* Left Accumulate: */ 

for j ^ 1 to Amax ddo 
fx^d + j 

M <~ Generainove{B\A\+s-d-j,d, fx, [s, fx]) 
end for 
/* RightAccumulate: */ 

for j ^ 1 to S — Amin do 
fn^ j 

T * GeiieralMove(Bs_j, /„, d, [/„, Amax]) 
end for 
end if 
return T 
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pegs used by disks from Bi during this process form an interval of contiguous integers, 
contained in the set of pegs available to this end, namely, the inclusion-wise maximal 
interval of pegs not occupied by any of the blocks Bi, . . . , -Bj_i. 

To move a block between pegs efficiently, all available pegs should usually be in use. 
More specifically, during the process of moving a sufficiently large block Bi, all of the 
available pegs are used. Furthermore, the algorithm allocates precisely h — i + 1 pegs to 
this end. This suggests that, in order to perform efficiently, the sizes of the h—1 blocks 
should satisfy rii > . . . > Uh-i = 1 (assuming n is sufficiently large). 

5.3.1 The Partition procedure 

In this section we present Partition — the procedure for partitioning a block B into the 
h — 1 blocks (Bi, B2, . . . , Bfi_i). We start by presenting an auxiliary function Remainder, 
which, for each stage j, provides the total number of disks to be assigned to the latter 
blocks - (-Bj+i, Bj^2, • • • , Bh^i). The definition of this function is given in Algorithm IH 

Algorithm 4 Remainder (ft,, n) 

/* Determines the size of the next set, by calculating the number of 'larger' disks. */ 
/* It is assumed that /i > 3 and n > 1. */ 
if n < ft then 

return max{n —1,1} 
else 

if ft = 4 then 

return min{n, round(\/2n)} 



Lemma 5.1 For any integers h > 3 and n > 1, we have 1 < Remainder(/i, n) < n. 
Furthermore, 

• lfh = 3, then Remainder(/i, n) = 1. 

• If h > A and 1 < n <h, then Remainder (/i, = n — 1. 

• If h > A and n > h, then Remainder (/i, n) >2. 
The proof is straightforward. 



else 



return min |n, 

end if 
end if 
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Algorithm 5 Partition(ft,, _B) 



/* Returns a partition of the n disks into h — 1 blocks of consecutive disks. */ 
/* It is assumed that h> 2 and i? is a non-empty block of disks. */ 
for j 1 to /i — 2 do 

ruj -s— Remainder(/i — j + 1, rij) 
Uj -S— rij — rrij 

Bj ^ [S,ni„, Bnun + Uj - 1] /* if n, = 0, then B, = */ 
B *^B-Bj 
end for 

/* B is now a singleton, so that Bh-i ^ {Bmax} */ 
Bh-i ^ B 

return (Bi,.. .,Bh-i) 



The formal description of the procedure Partition is given in Algorithm [51 
We argue that Partition is well-defined. To prove this, it suffices to show that at 
each of the h — 2 invocations of Remainder(/;, — j + 1 < j < /i — 2, we have 

h — j + 1 > 3 and nj > 1. The first of these inequalities follows from the fact that 
j < h — 2. Now observe that rii = n > 1, and n^+i = nij = Remainder(/i — j + 1, nj). 
Hence, by Lemma [STTl a simple inductive argument yields 

nj>l, l<j<h-2, (6) 

and we are done. 

In the following lemma, whose proof is straightforward, we collect for later reference 
a few properties of the partition (Bi, . . . , Bh^i). 

Lemma 5.2 The tuple (Bi, . . . , Bh_i) is a partition of B into blocl<s, satisfying: 

• -B/i-l = {-Bmax}- 

• \Bj\ = rij < n — 1 for each j E [h — 2]. 

• Each non-empty blocl< is lighter than all subsequent non-empty blocks in the parti- 
tion. 

It is easy to verify that, for a pair of indices i E [h — 1] and j E [h — i], 

Bj+,.i{h,B) = Bjih-i + l,Bii)), 

or, equivalently: 
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Lemma 5.3 For any integers 1 < i < h — 1 and 1 < j < h — i: 

nj+i_i{h, n) = nj{h — i + 1, n{i)). 

5.4 FarthestMove versus GeneralMove 

We assume without loss of generality that A = [l,h\. For any integers h > 3,n > 0, s 
and d, such that 1 < s < d < h, we denote by Gg^dih, n) the number of moves required 
by GeneralMove to move a block of size n from peg s to peg d using A. Similarly, we 
denote by F{h, n) the number of moves required by FarthestMove to move such a block 
from peg 1 to peg h. (Note that F = Gi^h-) 

It is easy to verify that, for h = 3, the algorithm GeneralMove works exactly as 
does the classical algorithm of [2Zl- In particular, it requires 3" — 1 moves to transfer n 
disks between the two farthest pegs in Paths, and ^^^^ moves to transfer them between 
neighboring pegs, yielding: 

Lemma 5.4 For any non-negative integer n: 

Gi^2(3,n) = ^F(3,n). 

5.4.1 Initial steps in the analysis of GeneralMove 

In this section we analyze the algorithm GeneralMove for moving a block B in Path/j, 
h > 3, from peg s to peg d, s < d, using the set A = [l,h] of available pegs. Let 
h' = s + h - d. (Note that h' < h.) 

Consider an index j G [s — 1]. At phase LeftSpread, a left-move of block Bj from 
peg s to peg j using h — j + 1 available pegs is performed, requiring Gj^s{h — j + 1, 
moves. Similarly, at phase Right Accumulate, a right-move of block Bj from peg 
j to peg d using h — j + 1 available pegs is performed, requiring Gj^d{h — j + l,nj) 
moves. 

Consider now an index j E [s,h' — 1]. At phase RightSpread, a far-move of 
block Bj from peg s to peg s + h — j using h — j + 1 available pegs is performed, 
requiring F{h — j + l,nj) moves. At phase LeftAccumulate, a left-move of block 
Bj from peg s + h — j to peg d using h — j + 1 available pegs is performed, requiring 
Gd^s+h-j{h -j + l,nj) moves. 

The remainder B{h') of blocks is moved in phase MoveRemainder, using a far- 
move from peg s to peg d, which requires F(h — h' + l,n(h')) moves. 
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The discussion above implies 
Lemma 5.5 

s-l 

h'-i 

+ J2 [^(^ - J + 1' ^j) + Gd^s+h-Ah -] + 1, n^)] 

j=s 

+ F{h-h' + l,n{h')). 
5.4.2 Initial steps in the analysis of FarthestMove 

In this section we analyze the algorithm FarthestMove for moving a block B from peg 
1 to peg h in Path/j, /i > 3, using the set A — [1, h] of available pegs. 

First, observe that the last block namely disk -Bmax, performs h—1 moves. 

Consider an index j e [h — 2\. At each of the phases Spread, Reverse, and 
Accumulate, a fcir-move of block Bj with h — j free pegs is performed, requiring 
a total of ?>F{h — j -\- 1, ^j) moves. Also, j — 1 shifts of block Bj with h — j + 1 free 
pegs are performed at phase Reverse, requiring altogether {j — l)Gi-^2{h — j + 
moves. 

For 1 < i < h — 1, denote by F{h,n)\n(i) the number of moves of the n{i) largest 
disks in the course of performing the algorithm FarthestMove. The explanation in the 
preceding paragraph yields 

Lemma 5.6 For any integers h>2,n>l, and 1 < i < h — 1, 

• F{h, n)\n(i) = Yl - i + 1> ^j) + U - i)Gi-,2{h + Uj)] + h-l. 

j=i 

i-1 

• F{h, n) = J2 [^F{h - J + 1, nj) + {j - l)Gi^2{h -J + l, n,)] + F{h, n) 

For the subsequent lemmas we put m = /i — A:+l,forl<A;</i — 1. 
Lemma 5.7 For any integers h > 2, n > 1, and 1 < k < h — 1, 

h~l 

F{h,n)\n(k} = F{m,n{k)) + {k-l)J2Gi^2ih-j + l,nj). 

j=k 
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Proof: By Lemma [5^ in the particular case i = 1, 

m-2 

F{m,n{k)) = y~^[3F(m — j + l,n-,(m, n{k))) 

+ (i - l)Gi^2{m - J + 1, nj{m, n{k)))] + h-k. 
By Lemma 15751 we obtain 

m-2 

F{m,n{k)) = ^[3F{m - j + l,nj+k-i) 
i=i 

+ {j - l)Gi^2im - J + l,nj+k-i)] + h-k (7) 

h-2 

5^[3F(/i - J + l,n,) + (j - k)Gi^2{h - J + 1, n,)] + h - k. 

j=k 

Observe that 6*1^2(2, n/i_i) = 1. Thus, by Lemma 15751 the right-hand side of (JT]) 
reduces to 

h-l 

F{h,n)\n{k) -{k- 1) ^G'i^2(/i -i + 

j=k 

and we are done. | 

Lemma IST^ and Lemma [5771 imply 

Corollary 5.8 For any integers h > 2,n > 1, and 1 < k < h — 1, 

k-l 

F{h,n) = Y.\-^F{h-j + l,n,) + {j-l)G,^2{h-J + l,n,)] 
i=i 

h-i 

+ F{m, n{k)) + (A; - 1) ^ Gi^2{h -3 + 1, 

j=k 

5.4.3 Moving from one End to the Other is the most Costly 

The following statement shows that GeneralMove requires the maximal number of moves 
when the source and destination pegs are at the extreme ends of the set A. 

Proposition 5.9 For integers h > 3,n > 1, s,d, such that l<s<d<h and d—s+1 < h, 

Gs^d{h,n) < F{h,n). 
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Proof: Denote h" = h — h' + 1. The proof is by induction on n, for all values of h > 3. 
For n = 1, we have nj = for each I < j < h — 2. Hence by Lemma 15.51 

Gs^d{h,l) = F{h",l) = h"-l = d-s < h-l = F{h,l). 

We assume that the statement holds for less than n disks and all h> 3, and prove it for 
n disks and all h > 3. Observe that 

1 < s + h-d = h! < h-l. 

By Lemma \5?2\ for each 1 < j < /i' — 1, we have nj < n. Thus, by Lemma [575] and the 
induction hypothesis, 

h'-l 

Gs^d{h,n) < Y,2F{h-j + l,n,) + F{h\n{h')) 

i=i 
h'-l 

i=i 

Since h' < h — 1 and Uh-i = 1: 

h-l 

Y,Gi^2ih-j + l,n,) > Gi^2{2,nH-i) = 1. 
j=h' 



Thus by Corollary 15.81 

h'-l h-l 

F{h,n) > ^3F(/i-j + l,n,) + F(/i",n(/i')) + (/i'-l)5^G'i^2(/^-J + l,n,-; 
j=i j=h' 

h'-l 

> Gs^d{h,n). I 



5.5 Proof of Theorem IXD 
5.5.1 Auxiliary statements 

Lemma 5.10 For any integers h> A and n>l, 

G^^2{Kn)<'^-F{Kn)-l. 
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Proof: By Lemma 15. 51 

h-2 

Gi^2{h,n) = Y,[F{h-j + l,nj) + G2^h-j+i{h-j + l,nj)] + F{2,n{h 
Note that F(2, n{h - 1)) = 1. Thus by Proposition \5M 

h-2 

Gi^2{h,n) < ^2F(/i-j + l,n,) + l. 
Note that h — 1 >3. By Corollary 15.81 in the particular case k = h — 1, 

h-2 

F{h,n) > ^3F(/i-i + + F(2,n(/i- 1)) + (/i-2)Gi^2(2,n/,_ 
i=i 

h-2 

> '^3F{h-j + l,nj) + l + h-2 



h-2 



> J2^F{h~j + l,nj) + 3. 

3=1 



Altogether, 



2 



I 



Lemma 5.11 For any integers h> A and n > 2, 



h-i 2 



Proof: By Lemma 15. 6[ 



h-2 



F{h, n) = Y, [^F{h - J + 1, n,) + (j - l)Gi^2ih - j + 1, n,)] + h - 1. 

Note that - 1) > |. Hence, 

2 2 f^'^ \ 2 

\i=i / i=i 
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We claim that 

Gi^2(/i-J + l,^,) < + l<3<h-3. (9) 

To this end, note that by Lemma IS.lOl this clearly holds if > 1. Otherwise nj = 0, so 

Gi^2{h - J + l,n,) = = '^F{h - J + l,n,). 

Since n > 2 and rih-i = 1, there exists a number j with 1 < j < h — 2 such that 

Uj > 1. 

The analysis splits into two cases. 
Case 1: nh-2 > 1- 
By dSl), 

h— 3 2 

i=i i=i 
Note that F(3,n,i_2) > 2. Thus, by Lemma EH 

1 _ 2 _ 1 

G'i_^2(3,?i/i-2) = -F(3,n/i_2) < -F(3,?2/,_2) - -. 



Altogether, we have 



^G'i^2(/i-J + l,%) = 5Z^i^2(/i-J + l,%) + G'i^2(3,n/,_2) + l 



i=i j=i 

2 2 2 

< -^F(/i-j + l,n,) + -F(3,n,_2) + - (10) 

/i-i 

3 



|5^F(/i-j + l,n,). 



By dH]), the right-hand side of (ITPI) is no greater than |F(/i, n), as required. 
Case 2; rij > 1 for some 2 M//?ere 1 < z < /i — 3. 
By (E]) and Lemma El 



^ Gi^2(/i-J + l,ri,) < 7T 5^ F{h-j + l,n, 



j&[h-2]\i je[h-2]\i 
By Lemma 15301 

_ 2 

0x^2 + 1,^0 < -F(/i-i + l,ni) - 1. 
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Altogether, we have 

h-l 

J2Gi^2{h-j + l,nj) = Yl Gl^2{h-j + l,nJ) + Gl^2{h-i + l,n^) + l 

By (IH]), the right-hand side of (JTI|1 is strictly less than ^F{h,n), and we are done. | 
Lemma 5.12 For any integers h > 5 and n > h, 

F{h, n) < 3F{h, ui) + y^l^ -l,n-ni). 



Proof: By Corollary 15.81 in the particular case k = 2, 

h-l 

F{h,n) = 3F(/i,ni) + F{h - l,n -n^) + Y,Gi^2{h - j + 1 

i=2 



We have n > h > 5. Thus, /z — 1 > 4, and by Lemma ISTTl we have n — rii > 2. Applying 
Lemma 15.111 with h — 1 and n — ni instead of h and n, respectively, we get 

h-2 2 

y^Gi^2{h-3,n^{h-l,n-ni)) < F{h-l,n-n^). (12) 

By Lemma 15.31 in the particular case i = 2, for each I < j < h — 2, 

n) = nj{h — l,n — rii). 

Consequently, 

h-l h-2 2 

^G'i^2(/i-j + l,^j) = YGi^2{h- j,nj{h-l,n-ni)) < -F{h - l,n -ni), 
i=2 j=i 

which provides the required result. 

Lemma 5.13 For any integers h > 5 and 1 < n < h, 

< = (7,,-n"'^ ■3^'^-"^ < U{h,n), 



where 9h and ah are as in Theorem \3.1\ and Ch = ■ 
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Proof: First note that U{h,n) < U{h,n), so it remains to prove that F{h,n) < 
U{h,n). Since h > 5 and n < h, Lemmas l5.ll and \5?2\ imply that rij = 1 for each 
i e [n-l]{J {h -1}. Thus by Lemma [5^ 

h-2 

F{h,n) = Y,[mh~j + l,n,) + {j-l)G,^2{h-J + l,n,)] + h-l 
i=i 

n-1 

= J2 - J + 1, 1) + (j - l)Gi^2{h -J + 1,1)] + h-1 
i=i 

n-1 

= ^3{h-j) + {j-l) + h-l = n{3h-n)-2h. 

It is easy to verify that for /i > 5 and n < /i, 

n{3h-n) - 2h < 3n(/i-2). 
Note that OhW^-^ > 1 for /i > 3 and n> 1. For a; > 1, we have x < . Therefore, 

and consequently, 

3n{h - 2) < ^^^n"'^ ■ 3^""'^ . 

Altogether, 

F{h,n) < 3n(/i - 2) < ^^^n"'' ■ 3^'^-"'^ = t/(/i,n). 

I 

5.5.2 Conclusion of the Proof 

The proof is by double induction on h > 3 and n. 

For h = 3, the algorithm works exactly as does the algorithm of |2Z] for the 3-m-a-row 
graph. Therefore, the number F(3,n) of moves required by this algorithm for n disks is 
3" — 1. The substitution h = 3 \n the upper bound U{h,n) suggested by the proposition 
yields: 

U{3,n) = Chu"'^ ■3^'^'''^^ = l-ri°-3i-" = 3" > F{3,n). 

For h = A, the algorithm works exactly as does the algorithm FourMove of Section l4.2l for 
moving n disks between the two farthest pegs in Path4. Therefore, as shown in Section 
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14.21 F(4,n) is bounded above by 1.6^/n ■ 3^^. The substitution /i = 4 in the upper 
bound U{h,n) suggested by the proposition yields: 



For h> 5 and n < h, Lemma 15.131 implies that F{h,n) < U{h,n). 

We assume that for arbitrary fixed h > 5 and n > h, F{h',n') < U{h',n') holds for 
all {h',n') with either h' < h or both h' = h and n' < n, and prove it for {h,n). 

Let m = nil = Remainder(/i, n). Then rii = n — m. By Lemma 15.121 

F{h, n) < 3F{h, n - m) + —F{h - 1, m). 

9 

The analysis splits into two cases. 
Case 1: n < ^^^—^ — 



(h-2)! ■ 

In this case, we have 



n < ^^n"\ 



{h-3)\ 



(13) 



and so, 



m 



Remainder(/i, n) = min < n 



((/i-2)!)°'^ 



-n 



Oh 



n. 



It follows that F{h,n) < ^■F{h — By the induction hypothesis, 

F{h,n) < ^U{h-l,n) = ilc,_i-n°'-^-3'^''-i-"'^, 

y y 



(14) 



By (112]), we have 



{h-3)\ 

Substituting (ITSi) in (IT4l) . we obtain 



3/1-3 



nh-3 



nh-2 



Oh-n—^. (15) 



F(/i, n) < — C;,_i ■ n"'^-! ■ a^'^-"^ . 
9 



(16) 



It is easy to verify that ^Ch-i < Ch and ah-i < oih- Thus we find that: 



Case 2: n > 
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In this case, we have 

n > 

and so, 

m = Remainder (/i, n) =min <|n, 
By the induction hypothesis. 



((/i-2)!)"^ 
((/i-2)!)°'> 



n 



Oh 



n 



Oh 



<n. 



F{h,n) < 3Ch{n - m)°''3^'^-("-'")''"' + —Ch-i ■ m'"^-' ■ 3«'^-l■"^'^-^ (17) 

9 



Observe that 



^ n J \ n 



Similarly, we have 

Ohin-m^jf^ = OhUh^ ( 1 



m\ — 



n 



{h - 2)n9tl 



(18) 



--enn^-l (19) 



and 



lh-3 



0tl 



1 

h-3 



3/1-3 



+ 1 



- , h-3 



(/i-3)(/i-2)! 

~ (h-3){/i-2)' 

in the range n > h > 5. Hence for n > /i > 5, we have 



(/i-3)(/i-2) 



1 

h-3 



fe-4 

n ''-2 . 



It is easy to verify that 'd{h,n) = (-fe_3)(fe_2) ^ ''"^ monotone decreasing with h and n 
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1/3 



?9(/i,n) < i9(5,5) = 
Put c* = (^)^^^ Now 

Substituting ^ and ((20]) in ([17]), we obtain 



(20) 



F(/l, n) < 3C/,n"'^ (l - ^^h-uT^-I + . . g^.-n^^+c* 

= C,,n°''(l-^)3^"-"'^ + ^3'=*C/,_i-m°''-i-3^''-"'^ (21) 
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The second and third terms on the right-hand side of (1211) may be omitted since: 



n 



Thus we conclude that: 



-1 1 
h—'i I 1 /1-2 ■m.'^-'i' 



-1 



h_3 \ h-3 

'h-1 



0. 



eh-nh-2 
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